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ABSTRACT 

We identify an abundant population of extreme emission line galaxies (EELGs) at redshift z ~ 

1.7 in the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) imaging from 
Hubble Space Telescope/ Wide Field Camera 3 (HST/WFC3). 69 EELG candidates are selected by the 
large contribution of exceptionally bright emission lines to their near-infrared broad-band magnitudes. 
Supported by spectroscopic confirmation of strong [OIII] emission lines - with rest-frame equivalent 
widths ~ 1000 A - in the four candidates that have HST/WFC3 grism observations, we conclude 
that these objects are galaxies with ^ 10 8 M 0 in stellar mass, undergoing an enormous starburst 
phase with M*/M* of only ^ 15 Myr. These bursts may cause outflows that are strong enough to 
produce cored dark matter profiles in low-mass galaxies. The individual star formation rates and the 
co-moving number density (3.7 x 10 -4 Mpc -3 ) can produce in Gyr much of the stellar mass density 
that is presently contained in 10 8 — 10 9 M 0 dwarf galaxies. Therefore, our observations provide a 
strong indication that many or even most of the stars in present-day dwarf galaxies formed in strong, 
short-lived bursts, mostly at z > 1. 


1. INTRODUCTION 
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The formation history of dwarf galaxies with masses ^ 
10 8 Mq can usually only be studied through ’archaeolog- 
ical’ age reconstruction, based on resolved stellar popula- 
tions (e.g., Grebel 1997; Mateo 1998; Weisz et al. 2011). 
Their high-redshift progenitors have so far remained elu- 
sive despite the ever increasing depth of spectroscopic 
observing campaigns and imaging from the ground and 
the Hubble Space Telescope (HST). In this paper we iden- 
tify an abundant population of z > 1 dwarf galaxies un- 
dergoing extreme starbursts, through HST / Wide Field 
Camera 3 (WFC3) imaging from the Cosmic Assem- 
bly Near-IR Deep Extragalactic Legacy Survey (CAN- 
DELS, Grogin et al. 2011; Koekemoer et al. 2011), that 
may well be the progenitors of present-day dwarf galaxies 
with stellar masses 10 8 - 10 9 Mq. 

At the present day, starbursts contribute a minority 
to the total star formation activity in dwarf galaxies 
(Lee et al. 2009). However, there is abundant evidence 
that the star formation histories are complex and that 
bursts play an important role (as reviewed by Mateo 
1998). Many authors find evidence for short-lived (^10 
Myr) SF events in nearby star-forming dwarf galaxies 
from a range of observational and modeling techniques 
(e.g., Schaerer et al. 1999; Mas- Hesse & Kunth 1999; 
Thornley et al. 2000; Tremonti et al. 2001; Harris et al. 
2004), while others argue that star formation epochs of 
dwarf galaxies are more prolonged (e.g., Calzetti et al. 
1997; Lee 2008; McQuinn et al. 2009). Simulations also 
indicate that star formation histories of low-mass galax- 
ies are episodic or even burst-like (e.g., Pelupessy et al. 
2004; Stinson et al. 2007; Nagamine 2010). 

As most stars in dwarf galaxies formed more than 5 
Gyr ago (e.g., Dolphin et al. 2005; Weisz et al. 2011), it 
is crucial to understand the mode of star formation in 
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Figure 1. Observed I — J vs. J — H colors (in AB magnitudes) 
from HST/WFC3 and ACS imaging for all objects in the UDS 
and GSD (small black points) and the sample of emission-line 
dominated objects (large red points with error bars), selected by 
I - J > 0.44 + a(I - J) and J - H < -0.44 - a(J - if), where 
cr(/ — J) and a(J — if) are the la uncertainties on the colors. The 
blue line represents the redshifted (z = 1.7) continuum colors of 
the Starburst99 model (Leitherer et al. 1999) for continuous star 
formation, in the age range from 1 Myr to 100 Myr. The red 
line represents the same model, but with the J-band flux density 
increased by the emission line luminosity predicted by the model 
(Starburst99 predicts Ha luminosity - [OIII] emission, which falls 
in the J band at z = 1.7, is assumed to have the same equivalent 
width). The black arrow indicates dust attenuation. 

dwarf galaxies at those early epochs, but ’archaeologi- 
cal’ studies do not have the resolution in terms of stellar 
population age to constrain strengths, durations, and fre- 
quency of bursts. The increased frequency of interactions 
with other galaxies and higher gas fractions at z > 1 may 
have resulted in strong, short-lived starbursts. In this 
paper we place the first constraints on the open ques- 
tion of how many and how frequently strong, short-lived 
starbursts occur in dwarf galaxies at z > 1, and how rel- 
evant this mode of star formation is for the build-up of 
the dwarf galaxy population in a cosmological context. 

2. DATA 

2.1. Multi- Wavelength Imaging 

We select objects from multi- wavelength photometry 
of two fields with HST/WFC3 and Advanced Camera 
for Surveys (ACS) coverage: the Ultra Deep Survey 
(UDS) field and the GOODS-South Deep (GSD) field 
at 4-epoch depth. 20 For the UDS we use WFC3 imag- 
ing in F125W ( J) and F160W (H) and ACS imaging in 
F814W (/) from CANDELS. For the GSD we use the J 
and H band imaging from CANDELS, supplemented by 

20 GOODS is the Great Observatories Origins Deep Sur- 
vey. CANDELS provides deep images over the central parts 
of GOODS-North and GOODS-South, and wider, less deep 
imaging over the remainder of those fields and over the other 
CANDELS fields, including UDS. See the CANDELS website, 
http://candels.ucolick.org/ , for details of the field layouts. 


WFC3 imaging from the Early Release Science (ERS) 
program (Windhorst et al. 2011), and I band imaging 
from GOODS (Giavalisco et al. 2004). The total area 
with /-, J-, and iJ-band coverage used here is 279 square 
arcminutes. 

Sources are detected in the H band with SExtractor 
(Bertin & Arnouts 1996) and photometry is performed 
with TFIT (Laidler et al. 2007), which uses additional 
imaging data sets, ranging from U to 4.5/im to produce 
resolution-matched, multi- wavelength catalogs. The cat- 
alog construction is described in full by Guo et al. (in 
prep.). In addition, we use a version of GALAPAGOS 
(Haussler et al. 2007) adapted for CANDELS WFC3 
imaging to measure structural parameters (van der Wei 
et al., in prep.). 

2.2. Color- Color Selection 

We select objects that are red in I — J and blue in 
J — H (see Figure 1), tracing luminous emission lines 
that contribute significantly to the total J-band light. No 
known continuum emission can produce such broad-band 
colors. The highlighted objects in Figure 1 have I — J > 
0.44 + a (I — J ) and J — H < —0.44 — <j(J — iJ), where 
a refers to the color uncertainty; that is, we select those 
objects that are significantly more than 50% brighter in 
J than in both I and H. We identify 69 such objects, that 
is, there is 1 per ~4 square arcminutes. They range in 
magnitude from TJab = 24 to J7 ab = 27, with a median 
of H = 25.8 (see Table 1). We note that there is no 
gap in color-color space between the emission line galaxy 
candidates that we select and the general distribution; 
the selected objects are merely the most extreme outliers. 

In Figure 2 we show false-color composites of all 69 
candidates. These sources are typically compact, but 
not unresolved; their J- and iJ-band half-light radii from 
GALFIT are typically 0.1". A subset (~ 20%) are more 
extended or consist of multiple components. We show the 
U through 4.5 fim spectral energy distributions (SEDs) of 
a subset of the emission-line candidates in Figure 3. The 
SEDs are seen to be almost entirely flat in F v , or in terms 
of ultra-violet spectral slope they have /3 ~ — 2, where (3 
is defined as F\ = . The J band is a notable outlier 

from this SED shape for all these objects. 

3. EXTREMELY BRIGHT EMISSION LINES 
3.1. Photometric Constraints 

No known objects have continuum SEDs that resemble 
those shown in Figure 3; in particular, the extraordinar- 
ily blue J — H colors are difficult to explain by any radia- 
tive process. Our hypothesis is that the J-band excess 
is due to one or more emission lines. The implied equiv- 
alent widths in the observed frame are extraordinarily 
high: EW ~ 1500 - 3000A. 

Among the emission lines that can reach such extreme 
EWs, Lya and [Oil] are immediately ruled out because 
the implied high redshift would produce a strong break in 
the SEDs; the lack of such a break implies z < 2 for these 
objects. WFC3/UVIS observations (Windhorst et al. 
2011) provide UV photometry over the ERS area. The 
average color of those candidates in the ERS area is 
F275W — U — 1.44, which suggests that the Lyman 
break is situated at around 3000A in the observed frame, 
which, in in combination with the very blue continuum 
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Figure 2. False-color composites, created from HST /, J, and H band image cutouts of the 69 emission line galaxy candidates. The 
cutouts are 3" on a side, the pixel scale is 0.06 // , and the full- width half maximum resolution at the longest wavelength (the H band) is 
~ 0.18". The IDs correspond to those in Table 1. The sources are typically compact, although a subset of about 20% have more extended 
morphologies or feature multiple components. 


slopes redward of the U band, implies 2 > 1.5. Only one 
candidate with UVIS coverage has a F275W — U color 
consistent with that of a galaxy at z < 1.5. The implica- 
tion is that strong [OIII] emission at 4959A and 5007A 
provide the most plausible explanation for the J band 
excess light. 

If [OIII] is responsible for the J band excess the redshift 
upper limit is z = 1.8. Furthermore, because we select 
objects with blue J — H colors, the H cannot contain 
the bright Ha line, which implies z > 1.6. Thus, solely 
based on their photometric properties, we suggest that 
our candidates are strong [OIII] emitters in the redshift 
range 1.6 < z < 1.8. 

3.2. Spectroscopic Constraints 

The hypothesis that [OIII] emission at z ~ 1.7 ex- 
plains the J band excess light is strongly supported by 


spectroscopic observations. While none of the candi- 
dates have ground-based spectra, WFC3 grism observa- 
tions are available for small portions of the GSD (one 
pointing in the ERS field, Straughn et al. 2011) and the 
UDS (from the supernova follow-up program 12099, PI 
A. Riess). The available grism coverage overlaps with the 
positions of 4 candidates in our sample (1 in the ERS, 3 
in the UDS), and strong emission lines are detected in all 
4 cases. The spectra (Figure 4) all show bright emission 
lines in the wavelength range 1.3 — 1.4/im, whose com- 
bined fluxes are in agreement with the excess light seen 
in the J band. 

The lines in all 4 spectra are readily identified as [OIII] : 
the asymmetry of the bright line, always extended blue- 
ward, is due to the two components of the [OIII] line, 
at 5007A and at 4959A, where the latter is 3x fainter. 
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Figure 3. Broad-band SEDs of the 12 emission line galaxy can- 
didates selected from the GSD field. Units on the y-axis are ar- 
bitrary, and the SEDs are incrementally offset by 0.4 dex in the 
vertical direction for clarity, sorted by continuum slope, indicated 
by the dotted lines. The objects are characterized by flat SEDs 
in F u over the entire range from U band to H band. The J band 
noticeably deviates from this trend as the result of strong emission 
line contributions, mostly [OIII] at z ~ 1.7. The observed wave- 
lengths for 2 : = 1.7 of various emission lines are indicated by the 
vertical dashed lines. 

In all cases, H ft is also detected. The redshifts are all 
in the range z = 1.65 — 1.80, in excellent agreement 
with what we inferred solely from photometry. We con- 
clude that our sample of extreme emission line galaxies 
(EELGs) form the high-EW tail of the general population 
of emission line galaxies seen in ACS and WFC3 spec- 
troscopic grism observations (e.g., Straughn et al. 2008, 
2009; Atek et al. 2010; Straughn et al. 2011). 

In principle, Ha emitters at 0.9 < z < 1.1 should also 
be included by our selection technique. The spectroscopy 
and UV photometry indicate that those must be far fewer 
than [OIII] emitters. Whether this is due to selection ef- 
fects or evolution in the number density of such objects 
remains to be seen. Still, even though our working hy- 
pothesis is that all 69 candidates are [OIII] emitters at 
z ~ 1.7, we should keep in mind that some fraction of 
our 69 candidates are likely Ha emitters at 2 ^ 1. 

3.3. Emission Lines and Broad- Band Photometry 

We have shown that selecting objects which are much 
brighter in J than in I and H works as a rather clean 
method for finding strong [OIII] emitters at 1.6 < z < 
1.8. Emission line galaxies with such excesses in other 
bands also exist, but a systematic search is more compli- 
cated as at most redshift ranges, multiple lines (most no- 
tably [OIII] and Ha) affect multiple photometric bands. 
Therefore, we refrain from conducting such a systematic 
search here. 

The existence of such emission-line dominated galax- 
ies complicates the interpretation of SEDs, which is es- 
pecially relevant in the case of the search for and SED 
modeling of rare, high-redshift objects. Although con- 


Figure 4. WFC3 grism spectra of the four candidates with grism 
coverage. The IDs refer to those in Table 1. GSD18 is object 402 
from Straughn et al. (2011); the 3 objects in the UDS are from 
supernova follow-up grism observations (program ID 12099, PI 
A. Riess). The three vertical dashed lines show positions of H/3, 
[OIII], and H a for 2 : = 1.7. These spectra strongly suggest that the 
majority of the objects in our sample are [OIII] emitters at 2 : ~ 1.7. 

tamination by emission lines is often considered to be 
a factor (e.g., Labbe et al. 2010), the extremely bright 
lines we observe suggests that their effect may be un- 
derestimated. Ono et al. (2010) explicitly showed that 
red colors in Lya emitters and z = 7 Lyman break 
galaxies may indicate the presence of evolved stellar 
populations or strong nebular emission lines (also see 
Schaerer et al. 1999; Finkelstein et al. 2011). Steep UV 
continuum slopes, such as observed in our objects, should 
serve as a warning sign for contamination by nebular 
emission lines at longer wavelengths to the point that 
those can dominate the broad-band flux density. 

4. STARBURSTING DWARF GALAXIES AT Z = 1.7 
4.1. Star formation or A GN? 

Before turning to our preferred starburst interpreta- 
tion, let us first point out that nuclear activity is not a 
likely explanation for the bright emission lines in the vast 
majority of EELGs. None of the objects in the CDFS 
have significant detections in X-ray or at 24/im. The ob- 
jects are spatially resolved in both J and H, and, more- 
over, the J and H band sizes are consistent with each 
other. 

Moreover, it is highly unlikely that all 69 objects 
are dominated by line emission from active galactic 
nuclei. At least at the present day, low-mass, low- 
metallicity AGN are exceedingly rare (Izotov & Thuan 
2008), much rarer than starbursting dwarf galaxies 
(Izotov et al. 2011). The implied black hole masses for 
the objects in our sample, as inferred from their UV con- 
tinuum luminosities (Shen et al. 2008) are 10 6 M© at 
most, when assuming an Eddington accretion of unity. 
At these low masses, at least at the present day, secular 
processes drive nuclear activity; thus, an unknown accre- 
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Figure 5. Rest-frame V-band absolute magnitude vs. rest-frame 
equivalent width of the [OIII] emission line at 5007A as calculated 
from the broad-band photometry as explained in the text and as- 
suming that all emission-line galaxy candidates are at z = 1.7. 
They span a range in luminosity, M v = —17 to My = —20, and 
have EW[oiii], 5007 between 500A and 1200A. See also Table 1. 

tion mode or triggering mechanism for nuclear activity 
would have to be invoked in order to explain an extreme 
change in the relative numbers of AGN and starburst 
powered emission-line dominated objects. At these low 
masses, merging cannot account for this. The starburst 
hypothesis, on the other hand, places these objects in the 
realm of dwarf galaxies, and their physical and statisti- 
cal properties are consistent with the abundances and 
masses of dwarf galaxies as we will discuss below. 

Although nuclear activity cannot be ruled out entirely 
- and line-strength gradients in star forming z — 2 galax- 
ies suggest that weak AGN may contribute to some ex- 
tent (Trump et al. 2011) - we can safely assume that the 
observed emission lines are effectively dominated by star 
formation activity. 

4.2. Starburst Ages and Masses 

We interpret the observations in the context of the 
Starburst99 model (SB99, Leitherer et al. 1999), which 
includes predictions for how the EWs of Hydrogen re- 
combination line evolve over time. Therefore, our first 
task is to estimate H/3 line strengths from the data. We 
attribute the excess light in the J band, compared to the 
continuum light measured in the I and H bands, to com- 
bined effect of emission lines. Therefore, we can compute 
the combined equivalent width as follows: 

Ew = ( tj — a) 

where Wj = 2845A is the effective width of the J-filter 
response curve, z = 1.7 to correct the observed EW to the 
rest frame, and / is the flux density f v in the respective 
filters. 

The relative contributions of the various emission lines 


are constrained by fitting Gaussian components to the 3 
emission lines seen in the grism spectra shown in Figure 
4, keeping the ratio between the two [OIII] components 
fixed at 3. We only use the 3 UDS spectra as H/3 is only 
marginally detected in the GSD spectrum. The emission 
line ratios are remarkably similar for all 3 objects: H /? 
contributes 1 /8 to the combined line luminosity, suggest- 
ing a very low metallicity (see, e.g., Salzer et al. 2005; 
Amorfn et al. 2010, for comparisons). Because the flux 
is dominated by the [OIII] 5007 line and is therefore more 
directly related to our observations, we show the inferred 
[OIII] 5007 EWs in Figure 5 (also see Table 1). How- 
ever, we model the observations by fitting the inferred 
H f3 EWs to the SB99 predictions. These are assumed to 
be always l/8th of the combined EW. The unavoidable 
intrinsic scatter in this conversion is mimicked by prop- 
agating a generous factor of two in the uncertainties of 
the quantities we infer below. 

EWu /3 is a sensitive age indicator, as it is quickly re- 
duced once a stellar population gains in mass or the star 
formation activity diminishes. For a SB99 model with 
continuous star formation with a Chabrier (2003) IMF 
with a high- mass cut off at 100M© and metallicity 0.2Z© 
the H /3 EWs imply that the galaxies in our sample typ- 
ically have ages of 10 — 20 Myr (Figure 6). It seems 
unlikely that these bursts will be much longer than this 
given the energy put into the interstellar medium (see 
Sec. 4.3); follow-up grism observations will directly con- 
strain the number of older bursts. 

If we assume a single burst model instead, we infer 
ages 3 — 5 Myr; all formation histories with declining 
star formation rates produce ages that are bracketed by 
these two extremes. Continuous star formation seems 
more realistic than an instantaneous burst, a notion that 
becomes physically untenable at very young ages. Star 
formation cannot happen faster than the dynamical time 
scale of -10 Myr for these systems. Most relevant for our 
analysis is that the stellar mass estimates derived from 
the two different models are very similar (see below). 

We note that choosing a different metallicity does not 
significantly change our results - a low metallicity is re- 
alistic for these low-mass systems (Amorfn et al. 2010), 
and the [OIIIJ-H/3 ratio suggests the metallicity is indeed 
low. Choosing an even lower metallicity (0.05Z© instead 
of 0.2 Zq) results in an increase in both age and mass by 
— 0.2 dex. We do suffer from the usual, unknown uncer- 
tainty due to our lack of knowledge of the stellar IMF. A 
different slope or cutoff at the high-mass end changes the 
number of ionizing photons from high-mass stars per unit 
stellar mass, and the entirely unconstrained number of 
low- mass stars (< 1 M©) determines the overall normal- 
ization of the stellar mass. In general, we conclude that 
the ages of these galaxies are < 40 Myr, which includes 
the intrinsic range in age and the systematic uncertainty 
due to the unknown star formation history. 

In the following we use the results from the continuous 
star formation model, but using the instantaneous burst 
model, by virtue of the insensitivity of the mass estimates 
to the choice of star formation history, does not change 
our interpretation and conclusions. 

For a given age, the SB99 model predicts the (rest- 
frame) V-band mass-to-light ratio, such that we can di- 
rectly estimate the mass after deriving the V-band lumi- 
nosity from the observed H band magnitude. We cor- 
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Figure 6. Masses, ages and star formation rates for the 69 
emission-line dominated objects in our sample, derived with the 
SB99 model, assuming that all are at z = 1.7, adopting a con- 
tinuous SF model with 0.2 times solar metallicity and a Chabrier 
IMF. The SFRs indicated by the diagonal lines are simply obtained 
by dividing the mass (x-axis) by the age (y-axis). The galaxies 
in our sample typically have 10 8 Mq stellar masses with young 
ages (5 — 30 Myr), or, equivalently, extremely high specific SFRs 
(~ 5 X 10 -8 yr-\ or ~ 50 x t~^ bble ). 

rect the luminosity and the derived mass estimate for 
extinction by comparing the continuum slope derived 
from the ACS photometry at rest-frame 2500A, typi- 
cally P 2500 ~ —2 (see Table 1), with the SB99 model 
prediction (rather constant at $2500 ^ —2.6 for the ages 
of these bursts). If we adopt the Calzetti et al. (2000) 
extinction law for starbursting galaxies the typical ex- 
tinction is E(B — V) = 0.2. 

The median mass we infer is 8 x 10 7 M 0 (see Table 
1 and Figure 6). Mass estimates inferred from the in- 
stantaneous burst model are only slightly smaller, by 
less than 0.1 dex on average. Internal consistency lends 
our modeling approach strong credibility: given the in- 
ferred ages, masses, and extinction corrections, the SB99 
model predicts dust-attenuated rest-frame UV luminosi- 
ties that are consistent with the observed rest-frame UV 
luminosities - the latter are not used in our modeling 
procedure. Thus, the model successfully describes the 
observed rest-frame UV and optical continuum spectral 
energy distributions as well as the observed emission line 
luminosities. 

Full modeling of the spectral energy distributions that 
includes emission line contributions will be presented in 
forthcoming studies that will also include objects with 
less prominent emission lines. As a consistency check 
with the results presented above we already applied 
the method outlined by Finkelstein et al. (2011) to the 
galaxies in our sample. This method compares the ob- 
served photometry with model spectral energy distribu- 
tions that include contributions from nebular and stel- 
lar continuum radiation as well as all nebular emission 


lines. Free parameters in this modeling procedure in- 
clude redshift, stellar mass, extinction and metallicity. 
We find photometric redshifts that are consistent with 
1.6 < z < 1.8 for the vast majority of the sample. More- 
over, the inferred stellar masses and ages are very simi- 
lar, even though this model is based on a different stellar 
population synthesis model. 

The implied star- formation rates (see Figure 6) may 
lead one to expect that these galaxies have significant 
24/im detections. However, we have verified that none of 
the 29 candidates in the GSD, which has very deep 24/im 
from GOODS, is detected. A possible explanation for 
this is the presumably low metallicity, which would result 
is relatively small dust masses and, hence, low infrared 
luminosities. 

The observed IRAC flux densities at 3.6/im and 4.5/im 
from GOODS and SEDS 21 are in most cases - there are 
two exceptions - fully consistent with the expected flux 
densities for the bursts observed in the UV. In addition, 
the galaxies have the same sizes in the J and H bands, in- 
dicating that the spatial extent of the region from which 
the line emission originates roughly follows the stellar 
light. Hence, there is no evidence for underlying older 
stellar populations. However, we cannot rule out their 
existence: maximally old stellar populations have mass- 
to-light ratios that are up to ^ 50 times larger than those 
of the bursts, even in the near infrared. 

If we assume a past star formation rate that is con- 
stant after averaging over >100 Myr time scales we find 
upper limits for the mass in older stars that is ^ 5 x the 
burst mass. The implied total stellar mass upper limits 
are then < 5 x 10 8 M 0 . This caveat notwithstanding, 
we assume in the remainder of this paper that there is 
no significant population of older stars in these galaxies, 
and that the observed bursts account for the total stellar 
mass. However, the bottom line is that the total stellar 
masses of these objects are well below 10 9 M 0 , in the 
regime of dwarf galaxies. 

5. DISCUSSION 

5.1. Comparison with Other Samples 

Galaxies with similar properties have previously been 
identified through broad-band photometry at z < 0.4 in 
the Sloan Digital Sky Survey (Cardamone et al. 2009), 
and have been shown by Amorfn et al. (2010) and 
Izotov et al. (2011) to constitute the most strongly star- 
forming and most metal poor tail of the well-known class 
of blue compact dwarf galaxies (e.g., Sargent & Searle 
1970; Thuan & Martin 1981; Griffith et al. 2011), which 
have very low metallicities and extremely high, spatially 
concentrated star-formation activity (Guzman et al. 
1998; Overzier et al. 2008). 

Cowie et al. (2011) (also see Scarlata et al. 2009)) 
recently studied the Lya properties of high-EW Ha 
emitters, providing a direct connection between higher- 
redshift searches of Lya (e.g., Ouchi et al. 2008; Hu et al. 
2010), and find Lya EWs ranging from 20 A to 200A. 
Combining this with the findings of Nilsson et al. (2011), 
who show that Lya emitters at z ~ 2 are objects with a 
very wide range in properties, it is clear that from Lya 

21 GOODS and SEDS (Spitzer Extended Deep Survey) provide 
the deepest IRAC imaging ever obtained, an excellent probe of 
stellar populations at high redshift 
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Table 1 

Sample of Extreme Emission Line Galaxies 


EELG2011 

RA (J2000) 
(deg) 

DEC (J2000) 
(deg) 

H 

(AB) 

EW lom] , 5007 

(A) 

/^2500 

log (M) 
(M 0 ) 

GSD1 

53.167064 

-27.858936 

24.67 

± 

0.07 

459 ± 40 

-1.83 

± 

0.09 

8.57 

± 

0.17 

GSD2 

53.080345 

-27.850572 

25.63 

db 

0.07 

569 db 67 

-2.04 

db 

0.19 

7.96 

± 

0.19 

GSD3 

53.046020 

-27.837322 

25.88 

db 

0.11 

507 ± 75 

-2.02 

db 

0.18 

7.94 

db 

0.20 

GSD4 

53.105087 

-27.819974 

26.20 

± 

0.11 

769 ± 143 

-1.75 

± 

0.30 

7.72 

± 

0.22 

GSD5 

53.067508 

-27.773595 

25.05 

db 

0.08 

566 db 74 

-2.19 

db 

0.19 

8.14 

db 

0.19 

GSD6 

53.097499 

-27.763919 

24.99 

db 

0.05 

700 db 53 

-2.12 

db 

0.09 

8.10 

db 

0.17 

GSD7 

53.122127 

-27.759542 

25.44 

± 

0.12 

535 db 99 

-1.70 

± 

0.28 

8.20 

± 

0.22 

GSD8 

53.171936 

-27.759145 

24.26 

db 

0.04 

693 ± 47 

-1.76 

db 

0.11 

8.52 

db 

0.17 

GSD9 

53.078754 

-27.750288 

24.86 

db 

0.04 

468 db 32 

-1.99 

db 

0.10 

8.42 

db 

0.17 

GSD10 

53.063690 

-27.745853 

26.37 

± 

0.09 

759 ± 134 

-1.56 

db 

0.33 

7.72 

± 

0.22 

GSD11 

53.007499 

-27.741867 

25.97 

db 

0.09 

534 db 76 

-2.13 

db 

0.19 

7.84 

db 

0.20 

GSD12 

53.114612 

-27.721979 

25.80 

db 

0.12 

641 db 139 

-1.94 

± 

0.31 

7.87 

± 

0.23 

GSD13 

53.101516 

-27.720882 

24.77 

db 

0.03 

490 db 29 

-2.35 

db 

0.08 

8.29 

db 

0.16 

GSD14 

53.055908 

-27.718803 

26.00 

± 

0.08 

501 ± 65 

-1.90 

± 

0.20 

7.94 

± 

0.19 

GSD15 

53.149536 

-27.710285 

26.64 

± 

0.15 

820 db 288 

-2.20 

db 

0.57 

7.36 

± 

0.30 

GSD16 

53.147617 

-27.707088 

26.10 

db 

0.08 

582 db 80 

-2.08 

db 

0.22 

7.75 

db 

0.20 

GSD17 

53.064220 

-27.706523 

25.41 

db 

0.05 

465 db 43 

-2.20 

db 

0.13 

8.12 

db 

0.17 

GSD18 

53.071292 

-27.705802 

25.24 

± 

0.04 

861 ± 66 

-2.36 

± 

0.11 

7.85 

± 

0.17 

GSD19 

53.181976 

-27.705038 

25.71 

db 

0.06 

1002 db 245 

-2.18 

db 

0.41 

7.72 

± 

0.23 

GSD20 

53.140815 

-27.692390 

26.23 

db 

0.10 

496 db 82 

-1.83 

db 

0.26 

7.88 

db 

0.21 

GSD21 

53.100936 

-27.676704 

24.76 

± 

0.10 

935 ± 139 

-0.95 

± 

0.28 

8.54 

± 

0.21 

GSD22 

53.118450 

-27.819919 

26.76 

± 

0.13 

870 ± 198 

-1.88 

± 

0.36 

7.41 

± 

0.24 

GSD23 

53.077606 

-27.812795 

26.81 

db 

0.17 

1512 db 338 

-2.22 

db 

0.28 

7.30 

db 

0.25 

GSD24 

53.132972 

-27.740102 

27.77 

db 

0.30 

698 db 318 

-2.23 

db 

0.47 

6.95 

db 

0.38 

GSD25 

53.084388 

-27.727920 

27.29 

± 

0.15 

562 db 164 

-2.39 

± 

0.43 

7.18 

± 

0.27 

GSD26 

53.141502 

-27.724880 

26.65 

db 

0.10 

650 db 106 

-2.10 

db 

0.25 

7.47 

db 

0.21 

GSD27 

53.112579 

-27.707090 

26.91 

db 

0.13 

954 db 262 

-2.40 

db 

0.42 

7.15 

db 

0.26 

GSD28 

53.046119 

-27.705604 

27.13 

± 

0.16 

1009 ± 293 

-2.14 

± 

0.45 

7.16 

± 

0.28 

GSD29 

53.139953 

-27.675138 

27.79 

± 

0.21 

1314 db 557 

-2.12 

± 

0.74 

6.93 

± 

0.37 

UDS1 

34.275299 

-5.274496 

25.38 

db 

0.09 

576 db 90 

-1.33 

db 

0.18 

8.31 

db 

0.20 

UDS2 

34.440769 

-5.262566 

25.74 

± 

0.09 

1081 db 147 

-1.41 

db 

0.17 

7.99 

=b 

0.20 

UDS3 

34.482173 

-5.261399 

25.28 

± 

0.08 

507 ± 95 

-1.88 

db 

0.24 

8.23 

± 

0.20 

UDS4 

34.268657 

-5.260064 

25.44 

db 

0.10 

614 db 83 

-1.53 

db 

0.14 

8.18 

db 

0.19 

UDS5 

34.426483 

-5.255770 

25.69 

± 

0.11 

701 db 95 

-1.66 

± 

0.10 

7.99 

db 

0.20 

UDS6 

34.428569 

-5.255318 

25.10 

± 

0.07 

731 db 86 

-2.12 

± 

0.13 

8.04 

± 

0.18 

UDS7 

34.325676 

-5.251743 

24.32 

db 

0.04 

656 db 43 

-1.59 

db 

0.05 

8.58 

db 

0.17 

UDS8 

34.314014 

-5.251047 

26.44 

db 

0.17 

728 db 153 

-1.39 

db 

0.20 

7.77 

db 

0.24 

UDS9 

34.382587 

-5.244620 

25.94 

± 

0.09 

478 db 64 

-1.82 

± 

0.14 

8.03 

d= 

0.19 

UDS10 

34.263534 

-5.239433 

25.48 

± 

0.07 

541 ± 64 

-1.84 

± 

0.14 

8.13 

± 

0.18 

UDS11 

34.311279 

-5.238957 

26.36 

± 

0.10 

735 db 94 

-2.43 

db 

0.12 

7.42 

db 

0.19 

UDS12 

34.473888 

-5.234232 

24.15 

± 

0.03 

713 db 42 

-1.72 

db 

0.05 

8.57 

d= 

0.16 

UDS13 

34.318141 

-5.232299 

25.35 

± 

0.07 

716 ± 68 

-2.12 

± 

0.08 

7.95 

=b 

0.18 

UDS14 

34.481567 

-5.222499 

25.69 

db 

0.11 

602 db 96 

-2.30 

db 

0.15 

7.81 

db 

0.20 

UDS15 

34.371166 

-5.214803 

25.45 

db 

0.09 

843 db 111 

-1.35 

db 

0.15 

8.14 

db 

0.19 

UDS16 

34.482921 

-5.214187 

25.39 

± 

0.08 

662 db 87 

-1.80 

± 

0.15 

8.07 

± 

0.19 

UDS17 

34.247516 

-5.205330 

25.95 

± 

0.10 

469 db 63 

-2.22 

db 

0.13 

7.89 

db 

0.19 

UDS18 

34.315448 

-5.200902 

25.83 

db 

0.13 

739 db 136 

-1.87 

db 

0.17 

7.84 

db 

0.21 

UDS19 

34.298866 

-5.191800 

26.25 

± 

0.15 

543 db 106 

-1.94 

=b 

0.17 

7.79 

=b 

0.23 

UDS20 

34.232082 

-5.190388 

25.16 

± 

0.06 

648 ± 55 

-2.20 

± 

0.06 

8.03 

db 

0.17 

UDS21 

34.308940 

-5.190090 

25.15 

db 

0.07 

609 db 66 

-2.28 

db 

0.11 

8.03 

db 

0.18 

UDS22 

34.416740 

-5.180443 

26.34 

± 

0.19 

1070 ± 307 

-1.72 

± 

0.31 

7.64 

± 

0.28 

UDS23 

34.387023 

-5.177240 

25.59 

± 

0.07 

591 db 74 

-2.15 

db 

0.14 

7.91 

± 

0.18 

UDS24 

34.252845 

-5.176362 

26.02 

db 

0.13 

779 ± 140 

-1.40 

db 

0.19 

7.91 

db 

0.22 

UDS25 

34.402145 

-5.175352 

24.47 

db 

0.05 

507 ± 39 

-2.11 

db 

0.06 

8.48 

db 

0.17 

UDS26 

34.459190 

-5.174448 

25.34 

± 

0.08 

552 db 59 

-2.29 

± 

0.09 

8.01 

± 

0.18 

UDS27 

34.284194 

-5.164084 

26.40 

± 

0.17 

576 ± 139 

-1.05 

db 

0.30 

8.01 

± 

0.26 

UDS28 

34.500236 

-5.155595 

25.89 

db 

0.14 

519 db 100 

-1.77 

db 

0.21 

8.01 

db 

0.22 

UDS29 

34.263267 

-5.152174 

26.35 

± 

0.13 

1003 ± 150 

-1.90 

± 

0.11 

7.57 

± 

0.21 

UDS30 

34.477771 

-5.147521 

25.66 

± 

0.14 

533 ± 103 

-0.65 

± 

0.21 

8.50 

± 

0.22 

UDS31 

34.296325 

-5.144416 

25.95 

db 

0.12 

546 db 88 

-2.17 

db 

0.14 

7.81 

db 

0.20 

UDS32 

34.419242 

-5.142892 

25.83 

db 

0.12 

721 db 138 

-1.69 

db 

0.21 

7.91 

db 

0.22 

UDS33 

34.246810 

-5.139120 

26.06 

zb 

0.11 

553 db 85 

-1.54 

=b 

0.19 

7.99 

± 

0.20 

UDS34 

34.371933 

-5.137272 

24.83 

db 

0.05 

586 db 50 

-2.73 

db 

0.08 

8.01 

db 

0.17 

UDS35 

34.314693 

-5.133675 

26.42 

db 

0.21 

1125 db 315 

-1.20 

db 

0.28 

7.80 

db 

0.28 

UDS36 

34.261756 

-5.134672 

25.32 

± 

0.08 

658 ± 76 

-2.16 

± 

0.12 

7.97 

± 

0.18 

UDS37 

34.380569 

-5.268105 

26.54 

± 

0.22 

832 ± 249 

-1.34 

± 

0.28 

7.71 

± 

0.29 

UDS38 

34.441444 

-5.215963 

26.60 

db 

0.19 

912 db 208 

-2.14 

db 

0.14 

7.38 

db 

0.25 

UDS39 

34.334854 

-5.177163 

26.99 

± 

0.18 

594 db 145 

-2.39 

± 

0.24 

7.26 

± 

0.26 

UDS40 

34.438095 

-5.160070 

26.86 

± 

0.17 

677 ± 157 

-2.29 

± 

0.21 

7.30 

± 

0.25 


Note. — EELG2011: identification number prefixed by the respective field acronyms; RA/DEC: coordinates 
from the CANDELS catalogs; H : H - band AB magnitude from the CANDELS catalog; -EW[oiii], 5007 : rest- 
frame equivalent width inferred from the /, J, and H broad-band photometry (see text for details); ^ 2500 '■ F\ 
continuum slope at rest-frame 2500A inferred from a linear fit to the B, V, and I broad-band photometry; 
log (M): stellar mass inferred from Starburst99 (Leitherer et al. 1999), as described in the text. 
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Extreme Emission Line Galaxies 


emitters one cannot derive a complete description of star 
formation in low-mass galaxies. On the other hand, Ly<a 
emitters at higher redshifts (z > 3) appear to be young, 
with small stellar masses (e.g., Finkelstein et al. 2009), 
similar to the emission line galaxies studied here. 

Narrow-band surveys identified galaxies with strong 
[OIII] and Ra emission lines with EW ~ 100 — 1000A at 
redshifts z = 0.3 — 1 (e.g., Kakazu et al. 2007), demon- 
strated to be young and metal poor (Hu et al. 2009). 
Most notably, Atek et al. (2010) pointed out the exis- 
tence of a class of emission line galaxies at z ~ 1.5 with 
EW > 1000 A that would most likely be included in our 
sample as well. However, so far, their nature has not been 
described and their cosmological relevance in the context 
of galaxy formation has remained unclear. Therefore, let 
us now put these starbursting dwarf galaxies in a cosmo- 
logical context. 

5.2. Cosmological Context : Implications for the 
Formation of Dwarf Galaxies 

Our sample with redshifts 1.6 < z < 1.8 consists of 69 
low-mass ( 10 8 Mq), young ( ^ 0.5— 4x 10 7 yr), extreme 

starbursting, presumably metal-poor galaxies. Their co- 
moving number density 22 is 3.7 x 10 -4 Mpc -3 , two or- 
ders of magnitude higher than that of nearby galaxies 
with similar EWs (Cardamone et al. 2009). The individ- 
ual star formation rates and the number density com- 
bine into 1.7 x 10 -3 Mq yr -1 Mpc -3 . This is a 2% 
contribution to the total star-formation rate density at 
z ~ 1.7 occurring in galaxies that contribute perhaps 
~ 0.1% to the total stellar mass density at that epoch 
(e.g., Karim et al. 2011). 

Placing these detections in the context of the burst 
fraction among equally massive dwarf galaxies is diffi- 
cult, however, since we cannot currently constrain their 
number density; galaxies without starbursts have older 
luminosity weighted ages and, as a consequence, are at 
least 2 magnitudes fainter than the starbur sting galaxies. 
With marginal detections, even in the latest WFC3 data, 
and no obvious spectral features, their redshifts cannot 
easily be estimated. Model predictions differ strongly 
from observational measurements of the galaxy stellar 
mass function (Guo et al. 2011) and cannot be used to 
constrain the burst fraction among the population of low- 
mass galaxies. 

Nevertheless, we can still gauge the importance of the 
observed burst for the formation of low-mass galaxies by 
making the reasonable (and testable) assumption that 
the observed bursts occur with the same frequency at all 
epochs l<z<3.5-a period of ^4 Gyr during which 
the cosmic star formation history peaked and after which 
the number density of starbursts declines as mentioned 
above. The basic indication that such bursts are im- 
portant is that the number of stars produced in such 
bursts over a period of several Gyr is comparable to the 
number stars in present-day dwarf galaxies. Let us use 
this consideration to construct a simple toy model that 
relates the observations presented here to the mass func- 
tion of present-day low-mass galaxies. Guo et al. (2011) 
use the data from Baldry et al. (2008) for the galaxies 

22 the values for our two widely separate fields, UDS and GSD, 
differ by only 12% 


with masses down to 10 7 M 0 , and their mass function 
can be represented by a simple power law for all galaxies 
with masses < 10 10 M©, that is, well below the knee of 
the Schechter (1976) function: 

/ M* \ -0-37 

« m *>~ 0 - 043 (hFm;) • < 2 > 

in units of Mpc -3 dlog(M*) -1 . 

Let us then express the stellar mass of the present- 
day descendants of the observed starbursting galaxies at 
z ~ 1.7 in terms of the following star formation history: 

ti /r M \) U1C s t x Nb urs t / 0 \ 

M desc = 7 ^ (A) 

/burst 

where Mb urs t ^ 2x 15 Myrx5 M e yr -1 = 1.5 x 10 8 M© is 
the total stellar mass produced in a single starburst (the 
factor 2 is included to convert the observed burst age to 
its total duration), TVburst is the number of bursts that 
occurs in each galaxy over the ^4 Gyr period between 
z = 1 and z = 3.5, and /burst is the fraction of the 
total stellar mass that is produced in such bursts - the 
rest is assumed to form in smaller bursts and/or a more 
quiescent mode of star formation. 

Now we can predict the number of present-day de- 
scendants by dividing the number of observed bursts at 
z ~ 1.7 by the duty cycle, that is, the fraction of the 
time that a bursts is occurring over the period of time 
that they can occur (here assumed to be 1 < z < 3.5, or 
4 Gyr). This duty cycle is the ratio of the duration of a 
single burst (2 x 15 Myr) and this 4 Gyr period, multi- 
plied by the number of bursts per galaxy A/burst- Thus 
we have 


0(Afdesc) ^ 


Mburst 

Aburst 


4 Gyr 
30 Myr’ 


(4) 


where Mburst is the co-moving number den- 
sity of the bursts observed at z ~ 1.7 

(3.7 x 10 -4 Mpc -3 dlog(M*) -1 , where we use the 
observed span in stellar masses of an order of magnitude 
- see Figure 6 - to introduce this differential unit). 

According to Eqs. 2 and 4 the mass of the descendant, 
Afdesc? is uniquely determined by the fraction of stars 
formed in bursts, /burst- Equating </>(Md eS c) (Eq. 4) and 
c t>(M *) (Eq. 2), and substituting for A/burst (Eq. 3) we 
obtain, after rearranging terms: 


Afdesc 
10 8 Mq 


2 4 f -1 - 6 

burst’ 


(5) 


Thus, if /burst is close to unity, that is, almost all stars 
are formed in bursts, then we infer that each galaxy must 
undergo one or two bursts on average and that, therefore, 
Mdesc ^ 1 — 2Mburst- It is perhaps more realistic to adopt 
a smaller value for /burst- If /burst ^ 0-5, we find that 
two or three bursts must occur in each galaxy, producing 
descendants with masses 10 9 Mq. The latter would 
imply a growth in stellar mass between ^ ^ 1.7 and the 
present by at least a factor 3 given the mass constraints 
on the underlying populations of the observed starburst 
galaxies. Galaxy formation models suggest that the typ- 
ical growth is indeed a factor of 3 or 4, but we should 
bear in mind that these models do not reproduce the ob- 
served co-moving number density evolution of low-mass 
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galaxies with redshift (e.g., Guo et al. 2011). Therefore, 
these predictions should be treated with care. 

Choosing /burst very low (< 0.1) implies a very 

large growth in mass, with high- mass descendants (> 
10 10 Mq). Such growth by, say, more than an order 
of magnitude in mass seems unlikely in the context of 
current galaxy formation models. In particular, models 
are better observationally constrained for these higher 
masses and the prediction is that more massive galaxies 
(10 10-11 Mq) grow in mass by a factor 3 or 4, not by two 
orders of magnitude, between z ^ 1.7 and the present 
(e.g., Keres et al. 2009). Given these constraints, the 
general conclusion we can draw is that our observations 
suggest that many or most stars in present-day dwarf 
galaxies (with masses < 10 9 Mq) have formed in a small 
number of starbursts at z > 1. 

The quantitative interpretation we offer here relies on 
our estimate of the characteristic burst duration (30 
Myr). However, the uncertainty in this assumed burst 
duration does not affect our conclusion that intense star- 
bursts at high redshift constitute an important phase in 
the mass buildup low-mass present-day galaxies. 

If the bursts last significantly longer, which we deem 
unlikely, then the total mass formed in a burst will be 
larger than assumed above, while the frequency of the 
bursts will remain the same since our selection method 
is only sensitive to young bursts. As a result, the fraction 
of stars formed in bursts, /burst will be larger than what 
we derived above. 

Significantly shorter bursts (say, ^5 — 10 Myr, as al- 
lowed by the SB99 model - see Sec. 4.2) are also un- 
likely, as the bursts are spatially extended, and the im- 
plied crossing time in these systems is about 30 Myr. 
Nonetheless, in the case that burst times are short, the 
burst frequency must be higher to account for their ob- 
served number (Eq. 4), but the total mass produced in 
each burst remains the same (see Sec. 4.2). As a result, 
the number of bursts per galaxy (Wmrst) and, therefore, 
the fraction of stars formed in bursts (/burst) increase. 

The main assumptions in our interpretation, then, are 
1) that the descendants of the observed galaxies remain 
low-mass galaxies up to the present day (< 10 9 M 0 in 
stars), and 2) that the observed bursts do not only oc- 
cur at z = 1.6 — 1.8, but are equally frequent over a 
much broader redshift range (1 < z < 3.5). The first as- 
sumption is supported by the understanding that more 
massive galaxies are not expected to grow in stellar mass 
by several orders of magnitude, as would be required if 
the descendants of the observed starbursting galaxies are 
much more massive. The second assumption is straight- 
forward to test observationally by searches for similar 
objects over a wider redshift range. 

5.3. Energy Budget and Core Formation 

The most remarkable property of these galaxies are 
their growth rates (specific star formation rates) of 20 — 
200 Gyr -1 . This is far outside the realm of normal, more 
massive star forming galaxies, which typically have ~ 
1 Gyr -1 (e.g., Karim et al. 2011). 

The amount of energy deposited into the interstellar 
medium through winds and super novae (10 56-57 erg) ex- 
ceeds the binding energy by an order of magnitude, even 
if the total mass is larger than the stellar mass by many 


factors. This implies that the gas that is fueling the 
starburst may be in the process of being blown out, re- 
sulting in the end of the starburst phase. Whether the 
gas will be blown out of the halo depends on the halo 
mass, which is currently unconstrained. If the gas does 
not leave the halo it may eventually once more cool and 
sink to the center of the potential well and trigger an- 
other starburst. Then, the observed bursts could be 
part of a semi-periodic cycle of star-formation activity, 
which is seen in simulations of low-mass galaxies (e.g., 
Stinson et al. 2007). Episodic star formation would also 
alleviate the difficulty of explaining the sudden occurence 
of such a large starburst in such small galaxies. If such 
starbursts are not cyclical, then another possible expla- 
nation, usually offered to explain the lack of star forma- 
tion at early times in even less massive systems, is that 
star formation had been suppressed at earlier epochs as 
a result of UV background radiation (e.g., Babul & Rees 
1992; Babul & Ferguson 1996). 

An intriguing possibility is that if most of the gas is 
expelled from the central region, then the stars them- 
selves could become unbound as well, dissolving the en- 
tire galaxy. Displacement of large amounts of mate- 
rial leads to changes in the potential, and if this hap- 
pens at very short time scales, then even the dark mat- 
ter profile can be altered. Outflows have long been ar- 
gued to play a role in producing cored density profiles 
in low-mass galaxies (e.g., Navarro et al. 1996). The 
time scales, masses, and star- formation rates we derived 
above for the galaxies in our sample roughly meet the 
requirements for such a process to be efficient. Recently, 
Pontzen & Governato (2011) showed that episodic star 
formation, even with less intense bursts, can produce the 
same effect. What the previous and subsequent evolution 
of the systems we observed is remains to be seen, but the 
energy balance and the possibility that such bursts are 
reoccurring make our observations consistent with this 
general picture; indeed outflows may be responsible for 
the formation of flattened dark matter profiles in low 
mass galaxies. These speculative scenarios can be tested 
further with observational constraints on the gas masses 
and hydro dynamical modeling of these systems. 

6. SUMMARY 

Our discovery of an abundant galaxy population at z ~ 
1.7 with extremely high emission line equivalent widths 
implies that many high-redshift, low-mass galaxies form 
many of their stars in extreme starbursts. We propose 
that we have observed an important formation mode for 
dwarf galaxies: a small number of strong starbursts that 
occur at early epochs (z > 1) each form ^ 10 8 M 0 in 
stars in a very short time span (~ 30 Myr) to build up 
the bulk of the stellar components of present-day dwarf 
galaxies. This is in quantitative agreement with ’archaeo- 
logical’ studies of present-day dwarf galaxies, which have 
shown that their star formation histories are burst-like 
and that the ages of their stellar populations suggest for- 
mation redshifts z > 1 (e.g., Weisz et al. 2011). Un- 
der the reasonable assumption based on ACDM predic- 
tions for galaxy grwoth that the observed galaxies grow 
in mass by less than an order of magnitude up to the 
present day, our observations provide direct evidence for 
such an early formation epoch and, in particular, that 
short-lived bursts contribute much or even the majority 
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of star formation in dwarf galaxies. 

Facilities: HST( ACS, WFC3) , VLT( VIMOS) , 

Spitzer(IRAC) 
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